A Dark Matter candidate from 
Lorentz Invar iance in 6D. 

Giacomo Cacciapaglia, Aldo Deandrea and Jeremie Llodra-Perez 

Universite de Lyon, F-69622 Lyon, France; Universite Lyon 1, Villeurbanne; 
CNRS/IN2P3, UMR5822, Institut de Physique Nucleaire de Lyon 
F-69622 Villeurbanne Cedex, France 
g . cacciapagliaOipnl . in2p3 . f r , 
deandreaOipnl . in2p3 . f r , 
j llodraOipnl . in2p3 . f r 

Abstract 

We study the unique 6 dimensional orbifold with chiral fermions where a stable dark 
matter candidate is present due to Lorentz invariance on the orbifold, with no addi- 
tional discrete symmetries imposed by hand. We propose a model of Universal Extra 
Dimensions where a scalar photon of few hundred GeV is a good candidate for dark 
matter. The spectrum of the model is characteristic of the geometry, and it has clear 
distinctive features compared to previous models of Kaluza-Klein dark matter. The 
5 dimensional limit of this model is the minimal model of natural Kaluza-Klein dark 
matter. Notwithstanding the low mass range preferred by cosmology, the model will 
be a challenge for the LHC due to the relatively small splitting between the states in 
the same KK level. 



1 Introduction 



Observations indicate that most of the matter in the universe is dark and of non-baryonic 
origin with a hint that it could be non-relativistic and weakly interacting at late times (cold 
dark matter). Various potential candidates have been suggested, among the most popular, 
there are weakly interacting massive particles (WIMPs). In the following we shall focus 
uniquely on this possibility. This group includes the lightest supersymmetric particle (LSP) 
obtained by imposing R-parity in supersymmetric models pQ; the lightest T-odd particle 
(LTP) obtained in little Higgs models with T-parity [2J; the lightest Kaluza-Klein particle 
(LKP) which is stable due to a residual parity which may survive in the effective theory 
if special conditions are required for the interactions at the fixed points [3j H] . Apart from 
phenomenological considerations, which render these symmetries practically compulsory in 
all these cases, the most satisfactory way to introduce a stable particle of this type is to 
consider a fundamental symmetry. From this point of view a residual Kaluza-Klein (KK) 
parity is the closest candidate, being a remnant of the extra dimensional Lorentz invariance 
after the breaking by compactification. However, the compactifications considered so far [3j[5] 
have fixed points or lines which correspond to lower dimensional objects (branes) where the 
extra Lorentz symmetry is absent at all. Therefore, the KK parity, which is automatic in 
the bulk of the extra dimensions, has to be imposed by hand on the branes. 

We consider these parities unsatisfactory and we want to insist on the idea of keeping 
a fundamental symmetry as the explanation of dark matter: is this possible? In the de- 
velopment of extra dimensional model building, the idea of branes and orbifolds with fixed 
points (both suggested by string theory) has proven itself very useful. For example, one 
can easily obtain chiral light fermions (from a theory which is inherently non-chiral in 4D). 
However, in order to preserve the KK parity, one needs to impose non trivial constraints on 
the fixed points. In models based on an interval, one requires the Lagrangian terms local- 
ized on the two physically independent end points to be identical: obviously, this is not a 
direct consequence of the compactification! Moreover, the presence of explicit breaking of 
Lorentz invariance introduces a large arbitrariness, often neglected, in the model, because 
the localized interaction terms only obey to a smaller subgroup of the extra symmetries. As 
an example, divergences which are forbidden in the bulk emerge again via localized counter- 
terms. The Ultra- Violet sensitivity of the theory is generically worsened opening a Pandora 
box of free parameters and potentially bad ultraviolet behaviour. On the other hand, the 
idea of large compact extra dimensions gives opportunities for model building allowing for 
elegant uses of fundamental symmetries, like in gauge-Higgs unification models [6], where 
the freedom to add localized terms is welcome in order to obtain a realistic model. 

For those reasons, we will require that Lorentz invariance is only broken globally by the 
compactification, and a KK parity is part of the orbifold, therefore exact and inevitable. In 
this spirit, the dark matter candidate is truly a remnant of Lorentz invariance, which we 
shall call by abuse of language, a lightest Lorentz particle (LLP) in the following. Moreover, 
all interactions are highly constrained, and the model will be relatively more predictive than 
standard KK parity models. The main hurdle to be passed is the requirement of a chiral 
spectrum for fermions, crucial to obtain the Standard Model in the low energy limit. 
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Orbifolds are quotient spaces of a manifold modulo a discrete group. The one-dimensional 
orbifolds are the circle S 1 = R/Z and the interval S' 1 /Z 2 : the circle has neither fixed points 
nor chiral fermions; the interval is the only orbifold with chiral massless fermions, however 
it possesses two fixed points (the boundaries of the interval). If we extend our four di- 
mensional world to include one extra dimension, therefore, there are no compactifications 
without boundaries which allow to obtain chiral fermions in 4D. The next step is to consider 
two extra dimensions. In the plane the possible isometries are translations (i), reflections, 
27r/n rotations (r) with n = 2,3,4,6 and glide-reflections (g), which are translations with a 
simultaneous mirror reflection. In 2D there are only 17 fundamental symmetry groups (see 
for example [7]) which correspond to the crystallographic groups in two dimensions, also 
called the wallpaper groups. The folding of the infinite periodic tiling of the plane can be 
described in a unique way by an orbifold. Only 3 of the resulting orbifolds are free of bound- 
ary or fixed points/lines. They are the torus IR 2 /pl, the Klein bottle M 2 /pg and the real 
projective plane M 2 /pgg. Our notation refers to the fact that space groups can be defined in 
a purely algebraic way: instead of specifying a representation of the generators, one can list 
the relations among them, so we have 

pi ~ Z 2 = (h) x (t 2 ) , 
Pg = (t 2 ,g\[g 2 ,t 2 }=0,t 2 gt 2 g~ 1 = 1}DZ\ 
Pgg = (r,g\r 2 = (g 2 r) 2 = l)DZ 2 ,Z 2 . 

Only the projective plane lR 2 /pgg allows chiral zero modes for fermions because it contains 
a 7r-rotation, as we will discuss in the following. We will therefore consider this geometry 
as the background for our model. One may push this exercise further and consider three 
or more extra dimensions. We will not do it in the following for the simple reason that 
the number of possible orbifolds increases very fast. Furthermore, adding more dimensions 
will lower the effective cutoff of the theory and therefore reduce the validity of the effective 
theory. The projective plane is the simplest possibility and the most reasonable one for 
building an effective theory which may have a wide validity range in energy. Here we will 
consider the minimal case where the Standard Model is embedded in this geometry: the only 
new parameter, neglecting higher order operators, is the size of the extra dimensions which 
determines the mass scale of the KK resonances. The geometry also contains two singular 
points where 4 dimensional terms can be added, and will be required as counter-terms for 
loop divergences: we will consider the most general such Lagrangians, with the spirit that 
they are small corrections (typically one loop level). This geometry was already proposed in 
the context of Grand Unification models [8], but no explicit example has been constructed. 

2 Chiral fermions without fixed points: the real pro- 
jective plane 

The minimal dimension of a fermion \l/ in 6 dimensions is 8 (contrary to the 4 components 
in 4 and 5 D) j|: the Clifford algebra contains 6 8x8 Gamma matrices T 1 . . . T 6 . Moreover, 
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one can define 

r 7 = r^r 3 ! 4 ! 5 ! 6 (2.1) 

which allows to define two 6D chiralities via the projectors 

p± = \ (1 ± r 7 ) • (2.2) 

Therefore, the minimal spinor representation of the Lorentz group are 4-component chiral 
fermions ty± = P±ty. Each one of those 6D-chiral fields contains two 4D Weyl fermions of 
opposite 4D-chirality. In order to be more concrete, we will use the following representation 
for the Gamma matrices 

and consequently 

where Pl and Pr are the projectors on the 4D chiralities. In this basis, the 6D-chiral fermions 
can be written as (using the Weyl representation for the 4D Gamma matrices) 

•-(&)■ •-(£)■ (2 - 5) 



with 



^± = ( X q ) , ^ = ( l ) , (2-6) 



where ip are Dirac spinors, x an d ?7 are Weyl spinors. 
The 6D action for a massless 6D-chiral fermion is 



s± = J dx 5 j dx 6 z -{*±r a d a y ± - (9 a *±) r a ^ ± } = 

+ ^ [^l±7s(<95 T id 6 )ipR± + iJ R ±-f 5 (d 5 ± i9 6 )^ L± + /i.e.] | ; (2.7) 



and a mass term 



St, 



1 dx5 J 


j dxQ 


m| 


I dx5 J 


j dx e 
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The only difference between the two 6D-chiralities is a different sign in front of the x§ 
derivative. This feature will be important when discussing the parity properties of such 
fields. 

The real projective plane can be defined by a 7r-rotation r and a glide g, whose action on 
the coordinates is 



g:{° ° 1 d ■ ( 2 - 9 ) 

Note that r 2 = (g 2 r) 2 = 1; for any field, the allowed parities are therefore 

p r = ±l; p g = ±l. (2.10) 
In terms of r and g one can define two translations: 
2 1 x 5 ~ x 5 + 27ri? 5 



t 5 =^: ^ ■ — t ( } 2. , ~„ . (2.11) 

No Scherk-Schwarz phases can be defined on this geometry: all the fields are periodic. Note 
also that we can define another glide 

/ f x 5 ~ -x 5 + 7lR 5 



g' = gr:{* & 1 " b , (2.12) 

under which the fields have parity p g > = p g p r . The two radii R5 and Rq are in principle 
different: here for simplicity we will fix R§ = R§ = R, and assume R = 1 in all the formulas, 
except introducing it back when discussing the phenomenology of the model. The size of the 
radius will determine the overall mass scale for the KK modes, itikk — 1/R- 

Let us know discuss the 4D chirality of the fermions. The glide g changes the sign of one 



coordinate only, x 6 : in order to keep the action invariant, from eq. 2.7, we see that the two 
6D-chiralities are exchanged. Therefore, the glide requires to start with a non-chiral theory 
in 6D. Under the glide, a generic fermion transforms as: 

*(g(x))= Pg r g *(x), r 9 = r 6 r 7 . (2.13) 

The two 6D chiralities are exchanged up to a sign, therefore for both parities a non-chiral 
4D massless mode is allowed. For this reason the Klein bottle, defined by a glide and a 
translation, does not allow for chiral fermions. 

Under the rotation r both extra coordinates change sign, therefore the two 4D chiralities 
must have opposite parity: a zero mode is allowed only for one of the two 4D chiralities. A 
generic fermion transforms as: 

#(r(x)) = Pr T r ^(x) , r r = ir 5 r 6 r 7 ; (2.14) 

with this definition, p r = +1 corresponds to a left-handed zero mode and p r = —1 to 



a right-handed one. Because of eq. 2.10 for each bulk fermion there is a massless chiral 
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Figure 1: Fundamental domain for a torus (left) and real projective plane (right). 

fermion. The real projective plane is therefore the unique 6D orbifold which allows for chiral 
fermions without fixed points. Note also that a bulk mass term is odd under the rotation. 
For completeness, the action of the second glide g' on a fermion is defined by: 

V(g'(x)) = PgPr r gr *(x) , T gr = iT 5 . (2.15) 

Note finally that the presence of a massless mode for each bulk fermion implies that su- 
persymmetry cannot be completely broken in this background: in fact, the supersymmetry 
generator, which is a 6D spinor, obeys the same properties, and therefore an unbroken N=l 
supersymmetry always survives. 

2.1 KK parities and singularities on the real projective plane 

The real projective plane is non-orient able and has no boundaries, however there are two 
points with conical singularities, where localized counter-terms can be added in general. 
Nevertheless, a KK parity is still preserved without any further assumption. Two symmetries 
can be used to define the orbifold, for example a rotation and a glide. The rotation has 4 
un-equivalent fixed points: (0,0), (it, it), (0, tt) and (tt, 0) (see figure [T]). The glide, however, 
transforms (0,0) — > (tt,tt) and (0,7r) — > (tt,0) (and viceversa); therefore, no fixed point is 
present globally, and the the eventual localized interactions on the corners of the fundamental 
square must be identified in pairs. The identification is not a consequence of an ad-hoc global 
symmetry acting on the UV completion of the model, as it happens in other orbifolds, like 
the chiral square [5], but it is part of the orbifold itself: the two identified points are indeed 
the same point in the 6D space. A crucial consequence is that a discrete KK parity is left 
unbroken: it can be identified with a translation by (tt, tt) (combined with r, this is equivalent 
to a 7r-rotation around the center of the square (7r/2,7r/2), thus it is exactly the same KK 
parity as in [5]): 

/ x 5 ~ x 5 + tt , . 

Pkk ■ < , (2.16) 

[ X 6 ~ X 6 + TT 



5 



Under this transformation, two identified singular points are transformed one into the other, 
and a generic KK mode with momentum (k, I) along the extra directions will acquire a phase 
(— l) k+l : therefore, (1, 0) and (0, 1) identify the lightest tiers of odd particles thus containing 
a stable Dark Matter candidate, while the (1, 1) modes can generically decay into zero modes 
via localized interactions. Note that the theory also possesses another KK parity 

Vkk ■ \ _ _ _ , (2-17) 

^ x§ ~ xq 

under which the (1, 1) states would be odd and therefore stable. If p ' KK is unbroken, the two 
tiers (1, 0) and (0, 1) would also pick different parity, and therefore contain two independent 
candidates. However, this symmetry requires that the localized interactions on the two 
points are the same: this is true for terms induced by loops of bulk interactions, however a 
generic UV completion of the model would violate such symmetry. Not being a fundamental 
symmetry of the orbifold, we will discard it and assume that it is generally broken. 



3 The Standard Model on the real projective plane 

Here we will consider a simple extension of the Standard Model on the real projective plane, 
and study the spectrum of the KK modes. More complicated constructions are in principle 
possible, and we reserve to study them in future publications. For now we will study the SM 
gauge group SU(3) c x SU(2) H /X U(l)y with a single Higgs scalar doublet, and a 6D fermion 
for each chiral SM fermion (doublets Q and L, singlets U, D, E and possibly the singlet 
neutrino N). The lowest order Lagrangian will be the SM one, extended to 6 dimensions. As 
we will shortly see, to each SM fields it corresponds a tower of massive resonances organized 
in tiers of modes, labeled by two integers (/, k) which correspond to the discretized momenta 
along the extra directions, and the field content of each tier will crucially depend on the 
parities of the fields under the orbifold projection. At leading order, all the states in each 
tier are degenerate with mass determined by the two integers 

Splittings within the modes in each tier can be generated by three mechanisms: the Higgs 
vacuum expectation value (VEV), bulk interaction loop corrections and higher order oper- 
ators localized on the singular points. Here we will focus on the simplest case, where, due 
to the flatness of the metric along the extra coordinates, the Higgs VEV is constant: an 
important consequence is that, due to the orthogonality of wave functions of KK modes 
from different tiers, the Higgs mechanism will not mix the tiers. Therefore, the KK expan- 
sion remains valid and the masses will be shifted, independently on the spin of the field, 
by the SM mass m , according to the formula: m 2 k = I 2 + k 2 + wIq. On the other hand, 
the loop corrections do generate level-mixings: being small effects, for the spectrum we will 
limit ourselves to the leading corrections, therefore diagonal terms only. Note also that the 
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loop induced terms will respect the full global symmetries of the space, and therefore, as 
an example, no splitting and/or mixing between the (1,0) and (0, 1) levels will be induced, 
neither decays of the (1, 1) modes into SM particles. Localized terms are generally required 
as counter-terms for the loop divergences: however, they must respect less symmetries than 
the bulk loops. The only unbroken symmetry will be the KK parity Pkk- In the following 
we will assume that localized terms are small as they are required at one loop level, and we 
will limit ourselves to a leading order expansion in them. 

In this section, we will first study the tree level spectrum for generic scalars, gauge fields 
and fermions on the real projective plane, thus identifying the possible field content of each 
tier, and then include the effect of the Higgs mechanism and lowest order localized terms on 
SM fields. One loops results for the lightest tier will also be presented because they play a 
crucial role when discussing the Dark Matter relic abundance. 

3.1 Scalars 

The action for a scalar field $ is (omitting the integral along the un-compact 4 dimensions) 

^scalar = J dx $ dx $ {d a &d a $ - M 2 $t$} , (3.2) 

a — 1, . . . 6; which leads to the equation of motion (EOM) 

(d 2 5 + d 2 + p 2 - M 2 ) $ = 0, (3.3) 

where p 2 = —d^d^. After Fourier transforming along the two extra coordinates, the field 
can be expanded in a sum of KK modes, whose wave functions satisfy the above equation 
with p 2 replaced by the mass square of the mode. The solutions of this equation are usual 
combinations of sines and cosines (with frequencies determined by the periodicity at 2%). 
The wave functions can be labelled with the parities under the rotation and glide: 





Pr 


Vg 


spectrum 


cos kx$ 


cos Ixq 


+ 


(_!)*+' 


k,l > 


sin kx$ 


sin Ixq 


+ 




k,l > 


sin kx 5 


cos Ixq 






k > 0,1 > 


cos kx 5 


sin Ixq 




(_l)fc+H-l 


k > 0,/ > 



The masses are given by the formula 

m 2 kl = M 2 + k 2 + l 2 . (3.4) 

The mass eigenstates can be labeled by their parity assignment {p r ,Pg) and KK number 
(k,l). Here is a full classification of the modes (with normalized wave functions): 
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(M) 


Pkk 


(++) 


(+-) 


(-+) 


(— ) 


(0,0) 


+ 


j_ 








(0,2/) 


+ 


-A- cos2/a; 6 






-i- sin2/x 6 


(0,2/-l) 






^ cos(2/ - l)x 6 


± sin(2/ - l)x 8 




(2k,0) 


+ 


-A- cos 2/crr 5 




-A- sin 2/crr 5 




(2k- 1,0) 






^ cos(2/c - l)x 5 




7b Sin ( 2/C _ 1 ) X 5 


(/c, /)k+l even 


+ 


- COS COS Ixq 


- sin sin Z:r 6 


- sin kx§ cos /a^ 


- cos kx§ sin /x6 


(k, l)k+i odd 




^ sin /ca;5 sin /x 6 


^ cos A;x5 cos 


^ cos A;a;5 sin /x 6 


^ sin kx 5 cos Za;6 



3.2 Gauge fields 

The action for an abelian gauge field is (also valid at quadratic level for non-abelian gauge 
symmetries) 

C 2w r 1 1 

%mge = J dx b dx & [ - -F al3 F a ? - — (d^ - £(d 5 A 5 + d 6 A 6 )) 2 ] , (3.5) 

where F a p = d a Ap — dpA a) and the {-gauge fixing term is added to eliminate the mixing 
between A^ and the extra polarizations A 5 and A 6 . The equation of motion for the vector 
component is 

- PFp, - -d^A^ + (d 2 5 + dl)A v = (p 2 + d 2 + dl)A u = , (3.6) 

which is the same as for a scalar field, and we have assumed that each KK mode satisfies 
the usual 4D equation in {-gauge: 

-d^F, u --^A,=p 2 A,. (3.7) 

Once the parities are assigned, the spectrum and wave functions will be the same as for the 
scalar field (with M = 0). 

The A^-Aq scalar sector is more complicated: in fact, the massive vector modes acquire 
their longitudinal polarization by eating a tower of scalar components provided by a combi- 
nation of A 5 and Aq, while another combination will give rise to a single tower of physical 
scalar states. Moreover, the parities of the scalar components are determined by the fact 
that they are part of a 6D vector, therefore if the vector component A^ has parities (p r ,P g ), 
the parities of A 5 and A Q components are respectively (—p r ,P g ) and (— p r , —p g ). In generic 
{-gauge, the EOMs are 

-dlA 5 + td 5 (d 5 A 5 + d 6 A 6 ) + d 6 (d 6 A 5 -d 5 A 6 ) = 0, (3.8) 
-dlA 6 + £d e (d 5 A 5 + d 6 A 6 ) + d 5 (d 5 A e - d 6 A 5 ) = 0. (3.9) 

Here we will focus on two gauge choices: the Unitary gauge, where all non-physical degrees 
of freedom are removed, and the Feynman-'t Hooft gauge, which is more useful for loop 
calculations. 
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Feynman-'t Hooft gauge 

In the Feynman-'t Hooft gauge £ = 1, the equations of motion for A 5 and A 6 decouple: 

(9% + d> - d*)A 5;G = , (3.10) 

therefore the two components can be treated as two independent scalar fields with proper 
parities. The wave functions and masses are the same as in the scalar case presented in 
detail in Section 13.11 



Unitary gauge 

In the unitary gauge £ —>■ oo, the combination 

d 5 A 5 + d 6 A 6 = 



(3.11) 



The two fields are not independent, and one can therefore expand both fields on the same 
tower of 4D scalars Afk,i)' 



i) 



(3.12) 



with c?505 + (9606 — 0. Using the latter relation in eq.s (3.8-3.9), the two wave functions 
respect the usual EOM of a scalar field: 

(p 2 + <9 5 2 + <9 6 2 )0 5/6 = O; (3.13) 

spectra and wave functions are again the same as in the scalar case, with the additional 



constraint from eq. |3.11| In the following tables, we list in detail the masses and normalized 
wave functions for the 4 possible parity assignments. In the (++) case the gauge symmetry 
is unbroken: 



(M) 

(0,0) 



Pkk 



A 



(++) 
r 

2tt 



-t=- cos 21xr 

V Z7T 



A 



(-+) 



.4 



(— ) 



(0,2/) 



=g sin(2l - l)x 6 



(0,21-1) 



-kr- cos 2kx$ 

V Z7T 



(2Ar,0) 



^ sin(2fc - l)x 5 



(2k- 1,0) 



(M) 



k+1 even 



- cos kxK cos Ixp 



- sin kx$ sin/x 6 



sin kx 5 cos lx 6 



cos kxc, smlxf 



(kj) 



k+1 odd 



nVk 2 +l 2 



cos kx$ sin/x6 



sin kx* cos Zxf 



If the gauge symmetry is broken by the glide(s), case (H — ), there is no zero mode in the 
spectrum: 
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(0,2/) 
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-A- sin 2/x 6 

V27T 




(0,2/-l) 





-4- cos(2/ - l)a: 6 






(2Jfe,0) 


+ 






-A- sin 2/crc 5 

V 2"7T 


(2k- 1,0) 




^ cos(2A; - l)z 5 






(A?, /)k+l even 


+ 


- sin /cx5 sill/a^ 


— 7= m cos A;x5 sin /x 6 

7TVK +t 


"ttv 7 ^^ 8111 C ° S ^ 6 


(fcj Ok+1 odd 




- cos kx 5 cos Ixq 


^ ;2 sin A;x 5 cos /x 6 


~7^^//fcp C ° S kX5 SmlX6 



If the gauge symmetry is broken by the rotation (and the glide g'), then there is a zero mode 
living in the A 5 component: 



(k, I) 


Pkk 


4~ +) 


4++> 


,(+-) 

^6 


(0,0) 
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l 




(0,2/) 


+ 




-A- COS 2/Xfi 

V 2"7T 




(0,2/-l) 




7b sin ( 2/ - ^ 






(2k,0) 


+ 


-A- sin2/cx5 

V Z7T 






(2k- 1,0) 








7^ cos(2A; - l)x 5 


(kj Z)k+1 even 


+ 


- sin kx$ cos /x6 


^VP+P cos fca;5 cos ^ 6 


— t= „ sin A;x5 sin Ixq 

7TVK +t 


(k, /)k+l odd 




- cos kx§ sin /x 6 


^^^2 sin fcar 5 sm/x 6 


^y^p c °s fca* cos lx 6 



Finally, the gauge symmetry can be broken by both rotation and glide g; in this case the 
zero modes resides in A 6 : 



(k,l) 


Pkk 


A* 


4 + -' 


,(++) 


(0,0) 


+ 






1 


(0,2/) 


+ 








(0,2/-l) 






vfc cos(2/ - l)x 6 




(2fc,0) 


+ 






-A- cos 2kx$ 

V 27T 


(2k- 1,0) 




^sin(2A;-l)x5 






(kj /)k+l even 


+ 


- cos sin /x6 


^^r^ sin fcx 5 smlx 6 


^VF+F cos /cx 5 cos /x 6 


(k, /)k+l odd 




- sin /cx 5 cos Ixq 







3.3 Fermions 

The action for a 6D Dirac fermion in Section [2] leads to the following EOMs for the 4 
components: 

2«^<9 M x± + (ds T id 6 )fj-s- = 0, (3.14) 
i*^fj ± -(d 5 ±id 6 ) X ± = 0; (3.15) 

as usual, we expand each component in a tower of 4D Dirac fermions fi >r (the subscript /, r 
indicate the 4D chirality) satisfying the usual Dirac EOMs. Those first order equations can 
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be decoupled [10], and each component satisfies the same quadratic equation as the scalar 
field in the previous section. The solutions are usual combinations of sin and cos, and the 
first order EOMs relate the coefficients of the two 4D components. The exact form of the 
solutions will be determined by the parity assignments of the fields. The rotation gives a 
different parity to the two 4D-chiralities, and a fermion with p r = + (p r = — ) will have 
a left-handed (right-handed) zero mode. On the other hand, the glide will relate the two 
6D chiralities, so that the four wave functions are not independent: the value of the parity 
under the glide does not play any role on the zero mode spectrum and, as we will see, the 
only requirement is that the SM doublets and singlets have the same glide parity in order to 
allow Yukawa couplings with the bulk Higgs. 

For a left-handed fermion, case (+±), the KK modes are given by: 



(k,l) 




x+ 




X- 


v+ 




V- 




(0,0) 

(o,0 
(M) 




1 




± 1 




-^sin/x 6 
— 77- sin kxz. 

Z7T 






-1)^ 
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27T 
1 

27T 


cos lx & 
cos kx 5 


±( 
±( 


-l) l ±coslx 6 
— l) fe 77- cos kxz. 

/ Z7T ^ 


i 

T( 


sin Ixq 
sin kx$ 



while for both k,l ^ 0, there are 2 degenerate solutions for each level which can be parame- 
terized as 



V 



(a cos kx 5 cos lx 6 + b sin kx 5 sin Ixq) fi 
±(— \) k+l (c sin kx 5 cos lx 6 — d cos kx 5 sin lx 6 
±(— l) k+l (acos/cxs cos Ixq — b sin kx^ sin/x 6 

(csin kxn cos Ixq + d cos kx$ sin Za; 6 ) f r 



'fr 

)fl 



\ 



(3.16) 



where we can use the EOMs and normalization condition to fix the coefficients 



a 
b 



k cos a—il sin a 



y/2-rr 
sin a 

V2tt 



d 



y / 2vrVfc 2 +« 2 
k sin a— il cos a 

v / 27r V / P+F 



The two orthogonal states can be obtained by choosing a 
an arbitrary mixing angle. 

For a right-handed fermion, case (— ±): 



(3.17) 



9 and a = tt/2 + 8, where 9 is 



(k,l) 






X- 




V+ 




»7- 


(0,0) 

(0,0 
(M) 




— 77- sin/xfi 

Z7T w 

7^ sin 






x / 2tt 
"V 2tt 




1 




± 1 


i(- 
i(- 


sin 

sin kx§ 


2^71 

7T- COS Ixr 

Z7T 

77- cos kx^ 


i( 


2X/27T 

-l)^cos/a; 6 
-l) fe 7^r cos kx 5 



and, for k,l ^ 0: 



( (asmkx 5 cosIxq + bcoskx 5 sin/x 6 ) fi \ 

±(— l) k+l (ccoskx$ cosIxq — dsinkx, 

±(— l) k+l (asinfcxs cos Ixq — b cos kx 
\ (c cos kx$ cos /x6 + d sin /cx 5 sin Ixq) f r J 



). 1 jl njtju 5 sin Zxg) 
cos fcrs sin/xe) 

sin /tc! y r 



(3.18) 
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with 



cos a k cos a+il sin a 

~ V2tt ~ V2ttv^+P (a i q\ 

J sin a r _ k sin a+il cos a yo.laj 

Note that we have fixed the normalization of the wave functions in such a way that the 
mass of each (k, I) KK level is real, i.e. m^n = \Jk 2 + I 2 . 



3.4 6D loop corrections 

The degeneracy of each KK level is removed at loop level: a complete calculation of the 
shifts is however beyond the scope of this paper. In the following we will focus on the modes 
(n,0) and (0,n) with n odd, in fact this case covers the lightest tiers, and the result will 
be important in the next section to determine the nature of the Dark Matter candidate and 
estimate its relic abundance. In general, the loop contributions (that we generically label II) 
can be divided in 4 pieces 

IT = n r + p g Ii G + p g p r Tl G , + p r Tl R : (3.20) 

the first term is the contribution one would get from the same fields on a torus and, after 
renormalization of the bulk kinetic terms, it leaves a finite contribution. The other three 
terms correspond to the two glides and rotation, in the sense that their sign depends on the 
parities (p r ,P g ) of the fields running in the loop. The contribution of the two glides is finite 
because the glides do not have any fixed points where a counter-term could be localized. 
On the other hand, the rotation does generate divergences which can be cut-off by counter- 
terms localized on the four points left fixed by the rotation, i.e. the two singular points of the 
orbifold. The singularities will be equally spread on the two points, because of the extended 
global symmetries of the bulk interactions. In the next section we will discuss the generic 
structure of the counter-terms: for now we will limit ourselves to cutting off the momentum 
integral in the loop and compute the coefficient on the log divergent term. Note also that 
bulk interactions respect both pkk and p' KK , therefore there will be no mixing between the 
states (n, 0) and (0,n): in the following we will compute the diagonal corrections, as the off 
diagonal ones do generate sub-leading corrections to the spectrum. 



3.4.1 Gauge bosons 

The tiers under study contain a gauge-scalar for each SM gauge bosons. For a generic gauge 
group, the mass receives corrections from the diagrams listed in Fig. |2j We performed the 
calculation in 3 different ways: we employed a novel method consisting in expanding in KK 
modes only along one direction and using the resummed 5D propagator [TTJ in the sum, and 
we checked the result, when possible, using the expansions in winding modes [12] and in 6D 
KK modes [13] . The methods are summarized in Appendix [Aj where we explicitly detail 
the three methods in the computation of the scalar loop "f" in the figure. The result of the 
calculation is summarized in the following table: 
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A* 



Figure 2: One- loop radiative corrections to the gauge scalar self-energy: gauge (a-c), ghost (d), fermion 
(e) and scalar (f-g) loops. 



Sm 2 gauge scalars 




X P 9 


XPgPr 


xp r 


a 


5T 6 


5 • 7C(3) 


3.(7C(3) + Si(n)) 


3n 2 n 2 L 


b 








-12S 2 (n) 





c 


-n 


-3 • 7C(3) 


-(7C(3) + S s (n)) 


5n 2 ir 2 L 


d 








-2B 2 (n) 





e 


-8T 6 











f 


T 6 


7C(3) 


(7C(3) + S!(n)) 


n 2 ix 2 L 


g 








-4S 2 (n) 






The contributions in the table must be multiplied by a normalization factor f ^3^5, where 
g = 2^ is the effective 4D gauge coupling and C{r) is a gauge group factor (defined as 
Tr(t"t£) = C(r)5 a ' b for a field in the representation r of a non-abelian group running in the 
loop, and the charge squared for a U(l) ). In the formula, T 6 is the typical sum appearing 
in the torus compactification [13] 



Ta 



1 

7T 



E 



0,0^(0,0) 



(k 2 + l 



2\2 



1.92 



(3.21) 



log 



is the log divergence associated with the rotation, and the ra-dependent contri- 



butions Si 23 are small corrections listed in the Appendix and coming from heavier modes 
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running in the loop. Summing over the SM fields, the corrections are: 

Sm% = [-79T 6 + 14C(3) + ttVL + B 1 - 4B 2 ] , (3.22) 

a 2 

5m 2 w = [-39T 6 + 70C(3) + 17ttVL + 7B 1 - 32B 2 - 2B 3 ] , (3.23) 

5m 2 G = [-36T 6 + 84C(3) + 2An 2 n 2 L + 9B 1 - 425 2 - 3B 3 ] . (3.24) 

Numerically, for n=l, the corrections to the mass 5m = |— are: 

5m B R = (-1.4 + 0.1L) • 10" 3 = -0.00094, (3.25) 

5m w R = (-0.4 + 5.8L) • 10" 3 = 0.026, (3.26) 

5m G R = (+0.2 + 28L) • 10" 3 = 0.13; (3.27) 



where we use a s {M z ) = 0.118, a(M z ) = 1/127, sin 2 9 W = 0.23 and A = 10 (L = 4.6). 
3.4.2 Fermions 

The corrections to the fermionic Lagrangian for a generic KK mode can be written in general 
as: 

5C = a L i)^ Vil P L i) + a R ^ Vil P R i> - bij^j . (3.28) 

The wave function renormalizations (in general different for the left-handed and the right- 
handed components) can be re-absorbed by a field renormalization, so that the shift in the 
mass (at leading order in the corrections) is: 

Sin? — b — m n aL ^ aR . (3.29) 

In the following we will list the contribution of the gauge and scalar (Higgs) loops to the 
three terms, in units of \ ^^i^ for the gauge loops (where ^(r) = (iV 2 — 1)/2N for a 
fundamental of SU(N), and the charge squared for a U(l) ) and 1 16 ^. 4fl for the Higgs (where 
y is the effective Yukawa coupling), to be multiplied by the parity of the bosonic field in the 
loop: 



fermions 




XP 9 


Xp g p r 


Xp r 


n 2 a,L gauge 





2 ■ 7C(3) 


2-(-7C(3) + F 1 (n) + F 2 (n)) 





n 2 a R gauge 





2 • 7C(3) 


-2F 2 (n) 





n 2 b gauge 





An ■ 7C(3) 


— 6nF 2 (n) 


4n 3 7r 2 L 


n 2 ai scalar 





7C(3) 


-7C(3) + F 1 (n) + F 2 (n) 


n 2 n 2 L 


n 2 a R scalar 





7C(3) 


-F 2 (n) 





n 2 b scalar 





2n ■ 7C(3) 


-nF 2 (n) 


n 3 n 2 L 
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The n-dependent terms are listed in Appendix |Bj For a generic fermion in the fundamental 
representation of SU(2) weak and SU(3) color and with hypercharge Yp\ 

6m * = M {< 21C(3) + inVL ~ Fl ~ ^ { Yh ' 2 + I 3 ' + §*0 + 

+i (21C(3) + nVL - fi - 2F 2 ) y 2 F \ . (3.30) 

Numerically, at level n — 1, for each SM fermion we find: 



StufR 


Q 


U 


D 


L 


E 


N 


light gen.s 


0.075 


0.067 


0.065 


0.012 


0.004 





third gen 


0.081 


0.072 


0.065 


0.012 


0.004 






where we used the same numerical inputs as for the bosons, and yt op = 1. 
3.5 Localized operators 

The divergences corresponding to the rotation require the presence of counter-terms localized 
on the two singular points. Due to the symmetries of the bulk interactions, the counter-terms 
should be equal. However, here we will take a more general approach, and add different terms 
on the two singularities: in this way we will break all additional global symmetries. In order 
to leave glide-invariance explicit, we define two localization operators: 

5 = i (5(x 5 )5(x e ) + S(x 5 - 7r)5(x 6 - tt)) , (3.31) 
5 n = - (8(x 5 )8(x 6 - 7r) + S(x 5 - tt)5{xq)) , (3.32) 

and label the two singular points with a subscript for the point (0,0) = (tt, tt), and n for 
(0,7r) = (vr,0). 

In general the localized interactions must respect only 4-dimensional Lorentz invariance, 
and the residual gauge invariance. For a scalar field like the Higgs many terms can be added 
including a mass term: here we will neglect this case because of the many free parameters 
and the low phenomeno logical interest of the Higgs resonances in this model. For a gauge 
field the situation is much simpler due to gauge invariance: in order to preserve the residual 
gauge invariance along the extra coordinate, we will write the localized interactions in terms 
of stress-energy tensor components [T4j : 

r — » ( Tli TP 2 T2i T? 2 _l_ r5i T? 2 _l_ T&i TP 2 _l_ r56i TP TPV\ . (Q QQ\ 

~ A~2 V — T ~ ~2~ 56 ~2 5/i ~2~ 6m ~Y b ^ b § ) ' [6.66) 

where % = 0, 7r, and the cutoff suppression compensates for the dimension of the 6D fields (2 
for a boson and 5/2 for a spinor). For a standard model gauge boson, with parities (+, +), 
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notice that both A 5j6 and d^A^ vanish on the singular points, therefore F 5)1 = F e/1 = and 
F 56 = d 5 A 6 - d 6 A 5 : 

b = %i (- - r f^ - d ^ 2 ) ■ ( 3 - 34 ) 

The first term corrects the kinetic term of the vector bosons, and it also introduces mixing 
between modes with different (k, I). Here we will assume those terms to be small, of the same 
order as the 1-loop corrections: this is a reasonable assumptions because they are in fact 
counter-terms required by divergences at 1-loop and their coefficient is suppressed by the 
cut-off of the model. Therefore, most off-diagonal terms will give higher order corrections to 
the masses: this is not the case, however, for tiers that are degenerate like (k,l) and (l,k). 
In fact, exchanging the two extra direction is a good symmetry of the real projective plane 
(this is true for R5 = Rq only, for different radii the degeneracy is removed at tree level). 
When expanding the localized terms in KK modes, the 2x2 blocks will have equal entries, 
while the loop contribution will be such that the diagonal entries are equal: from this, we see 
that the block can be diagonalized by the sum and difference of the two states^] Therefore, 
we define 

(k,l)±(l,k) 



(k,l) z 



V2 



with Z > k ,l,k 



0,1, 



■DO 



and parameterize the correction to the kinetic term as: 



2i 



Sa + 



4vr 2 A 2 



In this new basis (r±± = r w ± r^): 



(3.35) 



(3.36) 





(0,0) 


(0,20+ 


(21,21) 


(2k,2l) + 


{21 -1,21- 1) 


(2k -1,21- 1) + 


(0,0) 


n+ 


2r 1+ 


2r 1+ 


2V2r 1+ 


2r x _ 


2y/2n- 


(0,2Z')+ 


2r 1+ 


4r 1+ 


4n+ 


4v/2r 1+ 


4ri_ 


4v/2ri_ 


(2/', 21') 


2r 1+ 


4r 1+ 


4r i+ 


4^2r 1+ 


4ri_ 


4V2ri_ 


(2*',2Z')+ 


2^2r 1+ 


4^2ri+ 


4^1+ 


8r 1+ 


4V2n- 


8ri_ 


(2Z'-1,2Z'- 1) 


2ri_ 


4r x _ 


4r x _ 


4V2n_ 


4r 1+ 


4^2r 1+ 


(2Jfe / -l,2Z / -l)+ 


2y/2r X - 


4V2ri_ 


4 v / 2ri_ 


8ri_ 


4v/2r 1+ 


8r 1+ 



while (0, 20-, (2ife, 2/)_, (2/c - 1, 2/ - 1)_, (2/c, 21 - 1)± and (2/c - 1, 2Z)± (with Z > jfe) are not 
affected. The correction to the zero mode will renormalize the gauge coupling: 



9l 



4tt 2 1 + 



4tt 2 A 2 



96 

Art 2 



1 _ f 10 + r l7T 
4tt 2 A 2 



(3.37) 



lr The situation is more complicated when k and i are part of a Pitagorean triple such that k 2 + I 2 



or quartet with fc + I 



in this cases 3 or 4 states will be degenerate. However, this situation 



only happens for relatively large integers, the smallest ones being (0, 5) — (3, 4), (5, 5) — (1, 7), (1,8) — (4, 7), 
(2, 9) — (6, 7), (0, 10) — (6, 8) and so on. Therefore, the first case happens for states of mass 5y/2 ~ 7, which is 
too close to the cutoff and therefore phenomenologically not interesting: for this reason we will not explore 
the possibility of Pitagorean triples any further. 
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while the diagonal entries affect the masses: 



m tk,i) 



1 



4tt 2 A 2 



+ 



...). 



(3.38) 



Note also that many states are not affected by the localized terms: this means that the 
quantum corrections to such states are finite at 1-loop. 

A similar analysis can be performed for the scalar components of the gauge fields: note 
that no further mixing between the vector and scalars is induced, therefore the tree level bulk 
gauge fixing term is still appropriate. In this case, it is the (k, l) ev , modes to be unaffected 
(in Unitary gauge), and the correction is a mass term 



In the sum-difference basis, we obtain 



(3.39) 



Sij 


(0,2/-l) + 


(0,2/- 1)_ 


(2k, 21 - 1) + 


(2k, 21 - 1)_ 


(0,2Z'-1)+ 


4r 27r 








4^2r 27r 


(0,2J'-1)_ 





4r 20 


4\/2r2o 





(2k', 21' - 1)+ 





4^20 


8r 20 





(2k',2V 


4v^r 27r 








8r 27r 



This discussion can be easily extended for fermions where, due to the vanishing of one of 
the two 4D chiral components, no mass term can be added and only operators of dimension 
6 (like in the gauge boson case) are relevant. 



3.6 Electroweak symmetry breaking: the Higgs VEV 

Another source of mass is the Higgs VEV: here we will assume that the Higgs is (+, +) and 
its bulk potential contains a negative mass, so that a constant VEV is generated for the 6D 
Higgs scalar field. Due to the flatness of the VEV along the extra coordinates, it does not 
induce mixing between different KK tiers, therefore the KK expansion presented before can 
still be utilized to describe mass eigenstates. At the level (0, 0), therefore, we obtain precisely 
the SM spectrum and the absence of mixing with higher modes also means that the model 
does not suffer from dangerous tree level corrections to precision electroweak observables 
like the S and T parameters. For heavy tiers, the Higgs VEV induces a mixing between the 
weak neutral bosons, in a similar fashion as in the SM: the mixing angle is large because the 
states are degenerate at tree level, and loop effects will change the value with respect to the 
SM Weinberg angle. In the fermion sector, only the top will be significantly affected, and 
the two Dirac fermions corresponding to the left and right-handed SM top will be mixed. 



3.6.1 Gauge bosons: general analysis 

The Higgs VEV introduces new mixing between the vectors and some scalar components 
that must be cancelled by a suitable gauge fixing term. For an abelian gauge group, the 
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gauge fixing term in eq. |3.5| is replaced by 

r 1 

*-£-gauge " 



2£ 



{dpA* - £(d 5 A 5 + d 6 A 6 - gvMY 



when the Higgs is expanded 



H 



V2 



(v 6 + h + i(f>o) 



The EOMs of the vector part are modified simply by the addition of a mass term rriy 
For the scalars v4 5i6 and 4>q, the new EOMs are: 

2 - m v )A 5 + d 6 (d 6 A 5 - d 5 A 6 ) + m v d 5 (j)o + £d 5 (d 5 A 5 + d 6 A 6 - m v o ) 
(p 2 - my)A 6 + d 5 (d 5 A 6 - d 6 A 5 ) + m y <9 6 o + £d 6 (d 5 A 5 + d 6 A 6 - m V (j>o) 
[p 2 + dl + d 2 )<f) - m v {d 5 A 5 + d 6 A 6 ) + £m v {d 5 A 5 + 8 6 A 6 - m v <f) ) 

In the Feynman-'t Hooft gauge (£ = 1), the EOMs decouple: 

A 5 
A 6 



(P 2 



m 



l + dl + dl) 







(3.40) 



(3.41) 



9&1 



0, (3.42) 
0, (3.43) 
. (3.44) 



(3.45) 



and we have three independent towers with masses m 2 



b (i k) — I 2 + k 2 + m v , the only difference 
being the parities under rotation and glide. In the Unitary gauge, the condition 



m v cj) = d 5 A 5 + d 6 A 6 



(3.46) 



holds. Imposing this condition on the EOMs, we obtain the same decoupled equations as 
in the Feynman-'t Hooft gauge, however the fields are not independent anymore. Note that 
one solution of the constraint has the form: 

d 5 A 5 + d 6 A 6 = and O = ; (3.47) 

those solutions correspond to the physical scalars described in the previous sections, the effect 
of the Higgs VEV only appears in the extra mass contribution m v . The second independent 
combination of states which satisfies the condition is a new physical scalar, mainly consisting 
of Higgs component: if the gauge symmetry is unbroken by the orbifold 



(M) 


Pkk 


4 ++ > 




A—) 


(0,0) 


+ 








(0,2/) 


+ 


-= — , 2i = = cos 2lx & 

v / 2vr v /(20 2 +m2 r 




^ *pL -sin2/x 6 


(0,2/-l) 










(2Jfe,0) 


+ 


— — 2k cos 2kx 5 


^ _sin2fca;5 




(2k- 1,0) 










(&, Z)k+1 even 


+ 


1+^ cos kx§ cos Ixq 


^ sin kx 5 cos lx 6 


^ cos kx 5 sin lx 6 


(k, /)k+l odd 




JjJ'j smkx 5 smlx 6 


^ cos kx 5 sin lx 6 


l™ v ki sin kx 5 cos lx 6 



where N k)l = ^/(k 2 + l 2 )(k 2 + I 2 + m 2 v ). 
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3.6.2 Electroweak gauge bosons at the (1,0) and (0,1) tiers 

For scalar gauge bosons in the tiers (n, 0) and (0, n) with n odd, neglecting the localized 
operators, the correction to the neutral boson mass can be written as: 

( B n ).( 6 ™w + m \ -tzn6 w f w ).(W*\ 

v n ' \ - tan 6 w m 2 w 5m 2 B + tan 2 6 w m 2 w J \ B n J y 1 

Note that neither the Higgs VEV nor the loop corrections mix the two degenerate tiers. 
Analogously to the SM, this matrix can be diagonalized by 



Z n \ = f cos8 n sm6 n \ ( W* 
A n J I -sin6» n cos0 n J \ B n 



(3.49) 



with mass eigenstates 
n 2 1 



m A n ,z n = ^ + 2 + + 6m2w 



=F J (m| + 6m% - 5m 2 w ) 2 - Am 2 w (5m 2 B - bm 2 w ) ; (3.50) 



and mixing angle 



m 2 7 — m 2 A + m 2 7 — 2mL + Sm 2 R — Sm^ A/ 

tan 9 n = — ^ ^ ^ ^ ^ ^ . 3.51 

2m w m z sin V w 

Note that the mixing angle would be equal to the SM Weinberg angle if 5m 2 B = Sm^; also, 
due to the fact that the loop corrections grow with the KK mass scale, for large ttikk the 
mixing angle becomes smaller. A plot of the mixing angle as a function of tukk for the 
lightest tier is presented in Figure [3| while the splittings are plotted in Figure |4| The effect 
of the localized kinetic terms is easy to include: due to the fact that the bulk contributions 
(loops and Higgs VEV) are the same for the (n, 0) and (0, n) tier, they stay diagonal also in 
the basis (0, n)±. Therefore, the mass eigenvalues and mixing angle are given by the same 

2 B,W 2 B,W 

formulas as before, with 5m 2 B w Sm 2 B w + for (0, n) + and 5m 2 B w — > 8m 2 B w + 

for (0, 7i)_. 

For the charged weak boson, one obtains: 

Tl 2 

m w+ = + Sm w + m w ; (3.52) 

the contribution of the localized kinetic terms can be added in the same way as for the 
neutral bosons. 

Similar formulas apply for the vector states in the (n, 0) and (0, n) tiers, with n even, and 
to the (k, I) level: however, one needs to take into account the loop-induced mixing between 
the vector and the scalars, therefore the gauge fixing must be redefined at 1-loop. 
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Figure 3: Left panel: mixing angle sin6*i between the weak gauge scalars as a function of txikk- For zero 
KK mass we obtain the SM value, for large masses the mixing angle vanishes. Right panel: mixing angle in 
the top sector. 

3.6.3 Fermions 

The Yukawa couplings are only relevant for he top, therefore we will focus on this case, while 
the other fermions can be obtained by a simple generalization. The Yukawa couplings can 
be written as: 



SVukawa = -J dx 5 dx 6 Y 6 ^ Q H^ V + h.c. = (3.53) 

= -J dx 5 dx 6 Y 6 [i] Q+ Hxu- + Vq~ h Xu+ + XQ+Hf]u~ + XQ-Hfj u+ ) + h.c. 

This term can only be written if p r {Q) = —p r (U). Plugging in the wave functions, and 
denoting by qfj'^ and ufjf the left and right-handed components of the (k, I) KK mode, we 
obtain the corrections to the masses. For the zero modes: 

r Pg{Q) + Pg{U) Y 6 V 6 _ (0i o) (0,0) , , (<y fa \ 

£ Y ukawa(0,0) = 2 ^/f qi r ' 

it is clear that this term is non- vanishing only if the two bulk fermions have the same parity 
under the glide p g (Q) = P g (U) = p g , and 

m top =p g -j=-. (3.55) 

For (/, 0) and (0, 1) modes, we obtain 

^Yukawa(l,0)-(0,1) = -(-l)W>p (Ql U r ~ 9r«z) + h.C. . (3.56) 

For the (k, I) modes, the situation is more complicated due to the presence of 2 degenerate 
states: the mass term can be written in general as 

£ Yu kawa(k,i) = -(-l) fc+ ^to P (cosaQCOsac/ - sinaQsinac/) 

r tfM k ' l) -e' l) u?' l) )+h.c.-, (3.57) 
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Figure 4: Mass splitting between the different states in the tier 1 as a function of tukk'- in black the scalar 
photon (LLP), in blue the W and Z ', in solid red the gluon, in green the leptons, in magenta the light quarks, 
in dashed red the tops. 



where the two choices olq/u = 9q/u,tt/2 + Oq/u label the four independent states, and Oq/u 
are arbitrary parameters. If we chose 9q = —9u, two sets of states decouple so that there 
are no off-diagonal mass entries, and the mass matrices reduce to 

£Yu k awa( k ,l) = -(-lf+Wp {t^U™ - q^U^) + h.C. . (3.58) 

This is therefore a general expression valid for all modes. To find the mass eigenstates, we 
need to take into account the loop corrections to the Q and U masses. For the (1,0) and 
(0, 1) tiers: 

£ m ass = - ( ® 0, ) • ( * + )•(*)+ ^ ■ ( 3 - 59 ) 

V ' V m *°P R + ° m U J \ u r J 

The mass eigenvalues are 

m ti/2 = + m to P + 6m Q + \ 5m Q ±b) + Smu + ^5mu =p B^j , (3.60) 

with 



2 

1 Sm + Srnr 



B=\l(j l + """ w ^ — ) Wop- (3.61) 
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The eigenstates are given by 



til,: 
t<ll- 



cos a t ± sin a t 
=F sin a t cos a t 



Ql,r 
ULr 



and 



tana 4 = B 



1 5m Q + 5mu 
R 2 



(3.62) 



(3.63) 



4 Dark Matter and collider phenomenology 

The KK mass, and the Higgs mass, are the main free parameters of the model: calculating the 
relic Dark Matter abundance in this model, one can pin down the cosmologically interesting 
range for the KK mass. However, this is nothing but an estimate, because the result is 
very sensitive to the model of Cosmology and values of the cosmological parameters. In this 
work, we will assume the standard model of Cosmology and use the approximate formulas 
in Ref. [I]. A novel feature with respect to previous works in 6D [H |J5] is the smaller 
splitting between the states in the lightest tiers of resonances; therefore, we cannot neglect 
the co-annihilation with other particles species in the tier [16J. An average annihilation 
cross section can therefore be used to estimate the freeze out temperature and the relic 
density, assuming that all the particle species will decay into the LLP after freeze out. In 
our calculation, we will neglect electroweak symmetry breaking effects besides the mixing 
angle in the weak gauge boson sector which plays an important role in the calculation due to 
the relatively small mass splitting between the dark matter candidate and the heavier weak 
gauge resonances. We also included the main annihilation cross section between all the states 
in the lightest tier, and assume that the localized kinetic terms are negligible (therefore we 
only included the loop and Higgs contribution to the splittings). A more detailed study is 
left for a future publication. The result for the relic abundance as a function of the KK mass 
is summarized in figure [5] taking into account the presence of two degenerate tiers (1, 0) and 
(0,1), we find 

200 GeV < m KK < 300 GeV ; (4.1) 

and a limit ttikk < 400 GeV from the over-closure of the Universe. Note that the two 
degenerate tiers may be split in the case of asymmetric radii: when one radius is smaller 
than the other by more than few percent (in particular if the difference in mass is larger that 
the freeze out temperature, which is typically of order 4% of the KK mass) the heavier tier 
does not contribute significantly to the relic abundance and we obtain a range 

300 GeV < m KK < 400 GeV . (4.2) 

Note also that the mass range can be drastically modified by the localized kinetic terms: for 
instance, lowering the scalar gluon mass will increase the average cross section, and therefore 
push the preferred mass range to higher values. However, the presence of such largish 
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Figure 5: Relic abundance of the Dark Matter candidate as a function of tukk for degenerate (1,0) and 
(0, 1) levels (solid blue), and for the asymmetric case (dashed red). 

terms will also change the splittings of other levels, and modify the branching fractions into 
SM particles: we will not pursue this possibility any further at present. It is interesting 
to compare the preferred range with the results in Ref. [15], where the chiral square is 
considered. On the chiral square, the splittings are larger than in our case, therefore the 
lightest state is to a good approximation a purely U(l) hypercharge gauge scalar: because 
of this, the annihilation cross section is smaller, and the preferred mass range is lighter (for 
comparison, for light Higgs m# <C ttikk, they find ttikk < 200 GeV). In our case, the mixing 
angle is non negligible, and therefore the annihilation cross sections are much larger due to 
the SU(2) interactions. On the other hand, the co-annihilation with leptons (both singlet 
and doublet) dilutes the cross section due to the large number of leptonic degrees of freedom. 
Scalar gluons and quarks do not play a significant role, due to their larger mass. 

The inclusion of corrections from the electroweak symmetry breaking, which are impor- 
tant due to the lightness of the preferred mass, and resonant annihilation via the Higgs or 
the (2,0) and (0,2) tiers Q can change significantly the range: the inclusion of those effects, 
together with the localized kinetic terms, in a complete numerical calculation of the relic 
abundance is a subject under study [TT] . 

Another important issue is the compatibility with electroweak precision measurements [3] : 
in fact, following previous calculations, this low mass range may be excluded once loop 
corrections to precision observables are taken into account. However, a detailed calculation 
in this specific model has not been performed (and we leave it for a future publication); 
moreover, the low cutoff (naively 10 times the KK mass) means that bulk higher order 

1 We thank Hiroki Yamashita for pointing this out to us. 
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operators, similar to the SM ones, cannot be neglected. Therefore, the prediction power of 
the model is very limited in this sector: this issue is common to all other models of KK Dark 
Matter; a more detailed analysis is nevertheless required. 

In the following we will take the low mass range (rrixK < 400 GeV) as a ballpark to discuss 
the phenomenology of the model. The main feature of this model compared to previous ones 
is the relatively small mass splitting among the particles in each tier. The states in the 
first tier will chain decay into the LLP via tree level bulk interactions: the small energy 
available for the SM decay products makes their observation at the LHC challenging. Below 
400 GeV, the heaviest particles are the tops, about 70 GeV heavier than the LLP. Due to 
the large mass of the top, the main decay channel will be t^ ) -> bW™ -> bW*A^°\ 
where the virtual W* converts into a pair of quarks or leptons. The second heaviest particle 
is the scalar gluon, and it decays mainly G^ 1 ' ' — > qq^ 1 ' ^ — > qqA^ 1 ' '; the quarks from light 
families (from the up to the bottom) decay q( 1,0 ) — > qA^'°\ The final states from strongly 
interacting particles, therefore, will always contain 1, 2 or 3 jets and missing transverse 
energy, however the energy of each jet will be rather small, around 20 GeV, making their 
observation impossible at the LHC. The scalar W and Z will mainly decay to heavy leptons 
(or quarks for small masses): -> -> lvA^\ Z^°) -> llW) -> l~l+A^'°\ while 

leptons will decay directly to the LLP ft 1 ' ) — > lA^ 1 ' '. In this case the final state contains 
leptons, however their typical energy will be small again, less than 20 GeV, therefore they 
will likely escape detection at CMS and Atlas. In the table we summarize the main decays 
in the mass range 200GeV < tukk < 400 GeV (where MET stands for missing transverse 
energy) : 





m x - m LLP 
in GeV 


decay mode 


final state 
+ MET 


t d,o) 


70 


bWM) 


bjj 




blv 


G d,0) 


40-70 


qq (i,o) 


33 


g(l,0) 


20-40 


qA m 


J 


H/(!.o) 


20 


IvM), ulM) 


Iv 




20 


llM) 


11 


/d,0) 


< 5 




I 


AM) 










The second level accessible at the LHC is the (1, 1) [18], with mass V2m KK = 300 -=-450 
GeV. Besides chain decays similar to the ones for the lightest tier, the particles in this level 
can also decay to SM particles via localized interactions, only if such interactions break the 
extra accidental KK parity. The latter is the only decay mode for the lightest state in such 
level, therefore, if the level is stable, the same phenomenology as the lightest tier applies 
here, and the observation of such states will be very hard at the LHC. If the direct decays to 
SM are possible, one should easily observe resonances without missing energy. The states in 
this level can also be singly produced via the same small couplings that induce the decays: 
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as those are free parameters in the model, the single-production cross sections and branching 
ratios cannot be predicted and we will not comment on this level any further. 

Going up in mass, the next tiers are the (2,0) and (0,2). In principle, they can decay 
via bulk interactions in two states from the lightest tier, therefore they would go to invisible 
particles at the LHC. However, because of the fact that the mass is equal to twice the mass 
of the lightest tier at three level, the possibility that such decays are kinematically open 
really depends on the loop and Higgs induced splittings, and, in general, those decays will 
be suppressed by the small phase space. The loops (and kinetic terms) will also induce 
decays directly into SM particles and single productions, therefore, neglecting the localized 
terms, cross sections and branching fractions for this level can be calculated. The decays 
into SM particles will make this level easy to observe in final states without missing energy 
and with many clear resonances. 

The first observable missing energy signal will therefore arise at the next level, (2, 1) and 
(1,2), with mass y/5m KK = 450 4- 650 GeV. They are odd under the KK parity, and due to 
their mass they can only decay via a loop to a tier-1 particle plus a SM one. Therefore, the 
signature for such states will be missing energy plus one SM particle with hundreds of GeV 
of energy. Higher modes will repeat this pattern of decays. 

Another interesting feature of this model is the possibility of rare but spectacular signals 
with only one SM particle plus missing energy: this is true if the lightest tier is completely 
unobservable. For example, one can produce two tier-1 gluons, which radiate a hard gluon: 
this will generate a mono-jet plus missing energy signature [19] . Similarly, one may produce 
via loop coupling a (0, 1) state and a (2, 1) state, the latter decaying into (0, 1) plus a SM 
state, the only visible particle in the final state. More rare but spectacular signals may 
involve a single charged lepton, however those will be extremely rare due to the weak cross 
section but effectively background free due to the apparent violation of the electric charge. 

The model is being implemented in the FeynRules package [20] and it will be made 
publicly available, which will allow to interface the Lagrangian of the model to many Monte 
Carlo tools, and therefore study in detail the phenomenology simply sketched in this section. 

4.1 The "5D" limit 

The symmetries of the real projective plane allow to define two different radii along the 
two directions (this is not possible on the chiral square). Doing so, we would break the 
symmetry that exchanges the two directions, and remove the degeneracy between levels 
(n, k) and (k, n). In particular one can take a "5D" limit by sending one of the radii, say Rq, 
to zero: in this limit, the model collapses to a 5D model where the KK parity is imposed 
on the boundaries thanks to the extended symmetry of the 6D completion. However, the 
effective model still differs from the 5D case in a crucial way: the limit removes all the modes 
that carry momentum along the x§ direction, so that the structure of the tiers will closely 
resemble the 5D case, however the field content of each level will be very different. This 
can be rephrased saying that in the 5D limit, the 6D physics does not completely decouple: 
even if the modes that carry momentum along the xq direction decouple from the physical 
spectrum, the field content and spectrum of the remaining physical particles is in general 
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different from the one of a 5D model. On general grounds, this is the consequence of the 
fact that the starting point geometry is not factorized with respect to the extra dimensions, 
giving rise to a "geometrical" non-decoupling of the particle content and spectrum of the 
theory. 

Numerically the splittings can be very different from the 6D case. Nevertheless, we can 
estimate an indicative mass range from the asymmetric case in Figure [5j where only one of 
the lightest tiers is taken into account: the preferred range is 300-400 GeV. The range is 
lighter that in the usual UED 5D scenario (around 500 -j- 600 GeV [I]) because of the different 
spin of the LLP, and because of the co-annihilation with singlet and doublet leptons. The 
phenomenology of the model will be similar to the one of the symmetric 6D one, except for 
the absence of some of the levels (note that the model will appear as a "5D" one at the LHC 
as long as 1/Rq > few TeV). 

This example shows that even the minimal 5D model is not unique, and its structure 
and phenomenology may depend crucially on the presence of more dimensions which are 
too small to be observed. Our model, in the asymmetric limit, is therefore the minimal 5D 
model where the KK parity is naturally present, and the spin of the LLP can be considered 
a prediction. 

5 Conclusions 

We presented a class of models of Dark Matter in 6D where the presence of a Dark Matter 
candidate does not follow by an ad-hoc discrete symmetry, but it is a direct consequence of 
the topology of the compactification. The real projective plane is in fact the unique orbifold 
in 6D that allows both for chiral fermions and the absence of fixed points/lines, the latter 
ensuring the presence of an unbroken parity. Such parity, a relic of 6D Lorentz invariance, 
is exact even after including the effect of generic higher order operators localized at the two 
singular points of the compact space. Even though the model looks very similar to previous 
proposals in 6D (the chiral square), the topology of the compact space crucially affects the 
loop corrections to the spectrum and the structure of localized operators. We computed the 
one loop splitting of the lightest level, and found that the Dark Matter candidate is a scalar 
photon, with a mixing angle smaller that the Weinberg angle. The main difference with the 
chiral square is that the splittings in mass within the level tend to be smaller, therefore one 
needs to take into account co-annihilation with leptons and weak gauge bosons and the weak 
mixing angle cannot be neglected. An estimate of the relic abundance leads to a mass range 
200 < rriKK < 300 GeV (300 < rriKK < 400 GeV in the asymmetric case). In this range, 
the splittings are such that the observation of the particles in the lightest tier is virtually 
impossible at the LHC as the SM decay products are too soft. Higher levels are however easy 
to observe, as they may decay into SM particles without missing energy. The first missing 
energy signals may come from the level (2, 1) with mass y/hm^K- The small splittings also 
allow for rare but spectacular events with apparent charge non-conservation, like for instance 
mono-jet or mono-lepton plus missing energy, due to the fact that the decay products of the 
other quark or lepton are too soft to be detected. 
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The model also has a limit where one of the two dimensions is much smaller that the 
other, so that effectively we have one extra dimension with the KK parity imposed by the 
6D completion. Even though only one extra dimension is visible, the phenomenology of the 
model is very different from the usual 5D case, because particle content and mass splitting 
are very different. In particular, the KK odd levels do not have vector fields but scalars 
instead, and the Higgs resonances are missing. The DM candidate would be therefore a 
scalar instead of a vector massive photon. This is an example of a model in extra dimensions 
where the presence of dimensions at energies above the TeV scale can affect drastically the 
phenomenology. Moreover, the 5D model we present is the minimal model of KK Dark 
Matter in 5D where the KK parity arises naturally, and the prediction is the different spin 
of the lightest stable particle. 

Acknowledgments 

We thank Michele Papucci for useful discussions and comments. We also thank the Ecole de 
Physique in Les Houches where the preliminary results of this work were presented during 
the workshop "Physics at TeV colliders" , II session, and were we had useful discussions with 
Renaud Bruneliere, Christophe Grojean, Hiroki Yamashida, Jay Wacker. The research is 
supported in part by the ANR project SUSYPHENO (ANR-06-JCJC-0038). 

A Appendix: loop corrections to the masses 

The calculation of the loop corrections to the masses can be performed using three methods. 
First, we can use the expansion of the 6D propagator in winding modes [12]: in this way 
it is straightforward to renormalize the 6D kinetic terms, which corresponds to removing 
the contribution of the zero winding modes [13J. However, the calculation is challenging 
in general due to the presence of Bessel functions in the expansion. The second way is 
to expand in KK modes along one direction, and use the resummed 5D propagator [TT] 
along the other: the advantage is clear when computing corrections to (n, 0) modes, where 
conservation of momentum along the second extra direction simplifies the sum over the KK 
number in the propagators. Finally, one can use the usual KK expansion [13]: in this case 
a more sophisticated technique is required to renormalize the kinetic term. 

As an illustration of the three techniques, we will detail the explicit calculation of the 
scalar loop "f " from Figure [2] for the A e scalar modes (n, 0) with n odd. The results from 
the other loops has been calculated using at least two of those techniques. 

A.l 6D winding modes method 

A 6D scalar field satisfies the following equation of motion : 

(-d M 0" + 0g + 0g)0 = O. (A.l) 
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It is convenient to calculate the propagator in a mixed momentum representation in the un- 
compactified 4D and position space along the extra directions. The propagator is therefore 
the Green function of the following operator 

(p 2 + d 2 + dl) Gf(jp, x 5 - x' 5 , x 6 - x' 6 ) = i8(x 5 - x' 5 )S{x 6 - x' 6 ) . (A.2) 
The solution [12] , defining p = a/p 2 and ~y = (x^^xq), is: 

G 6 s D (p, ? - 70 = \h£\ p \f - 7'|) , (A.3) 

where Hq is the zero order Hankel function of first kind. The propagator on the real 



projective plane is given by [21] 

+Pr G% D (p, V - r(V') + il ) + p r p g G 6 s D {p, V - r * gtf') + Q )] (A.4) 

where f2 = (2imi, 27m 2 ) with (nl, n2) G Z 2 forces translation invariance ((ni,n 2 ) are the 
winding modes), p r and p 9 are the parities of the scalar field under rotation and glide and 
f(y') are the transformed of the point y' under the transformation /. 
The loop correction to the A e propagator is given by 

TI 66 = 2^(01^^ 

= iN (2vr) 2 f d'k [ Gf(k, y>, ~y)At°\q, ~y)A ( t°\q, (A.5) 



where A^'°\q,~y) = ^^sinnxs is the wave function of the external field, r] 66 = — 1 is a 
metric factor and g\ = (27t) 2 g 2 is the 6D gauge coupling. In the following, in order to simplify 
the notation, we will always omit the normalization factor iV = 2g ^g^ 7 " 36 ■ The correction 
ill 66 can be split into four terms whose signs depend on the parity of the scalar field under 
the symmetries of the space: 

n 66 = n T + Pg n G + Pr u R + PgPr n G , ; (A.6) 

Ut is the contribution we would obtain on a torus of same radii and it is finite after the 
kinetic term renormalization, the other three terms are generated by the symmetries of the 
orbifold and we do expect a log divergence arising in 11^ due to the fixed points of the 
rotation. 
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Torus 



The torus contribution is given by the first term in eq. (A. 4) plugged in eq. (A. 5): 



N 

7 



n 



-4tt 3 / dk E Yl k% K (2irk E ^nl + nl) (A.7) 

^° (nl,n2)£Z 2 

where Kq is the K-Bessel function of zero order and we have performed the Wick rotation 
to write the last integral in Euclidean space. The zero winding mode (n 1; n 2 ) = (0,0) 
contribution is UV divergent, however such divergence is the same we would get in the limit 
of un-compactified space, therefore it can be absorbed by a wave function renormalization 
of the 6D field. Removing the (0, 0) mode from the sum and integrating in k 

4 



o 



dk E k\ K (k E a) = (A. 



we obtain: 



n r =£ T. with T 6 = i £ ^^1.92. (A.9) 

(nl,n2)^(0,0) V 1 27 



Glides 



From the second term in the propagator eq. (A. 4) 



n G = ^ / d*k Jd-^Yl \n$\k IV - gilt) + Tt \f^P . (A.10) 

As the glide does not change sign to the x*, component, the Hankel function does not depend 
on £5 and its integral will lead to the normalization of wave functions. After Wick rotating 
and integrating in k as before we obtain 

n r 27r 1 

n G = — V / dx 6 ? , (A.ll) 

4?r (m,^ J° (K " V2) 2 + (* 6 A + n 2 - l/2ff 

where we numerically checked that 
N 



Ti' 2 



(nl,n2)GZ 2 



X) / dx5 7r , ^2 , ( / , 77^2 =7C(3). (A.12) 



((m- 1/2) 2 + (x 6 /tt + n 2 - 1/2) 2 ) 



One obtains a similar expression for the other glide n G ', now x G can be easily integrated 
out and we are left with 

tt JL ST Fa l-cos2m; 5 

" A *\n^'° ((X5/7F + ^-V2) 2 + (n 2 -l/2) 2 )2 

= j(7C(3) + S!(n)), 
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where we again numerically checked that the n-dependent term corresponds to the function 
in Section E31 



Rotation 

The last contribution is coming from the rotation part: after the Wick rotation 

AT poo r y 

Tl R = — 2tt / dk d~yJ2 k3K o( k \ 2 ^+ Q \) sin2nx 5- ( A - 14 ) 
n 

To extract the divergent part, we cut-off the 4D momentum A; at a scale A, and numerically 
integrated over the compact space and summed. One can therefore show that the integral 
is equal to 



N 9 9 , A 2 + n 2 

n R =- nVlog — . A. 15 

4 n z 



From the integral form, one can see that the divergences appear when \2~y + Q, \ = 0: those 
points are indeed the fixed points of the rotation, i.e. the corners of the fundamental square. 
In his notation, their geometrical origin is clear. 

A. 2 6D mixed propagator method 

Using the full 6D propagator is complicated by the fact that one needs to deal with Bessel 
functions and re-sum a double sum. On the other hand, 5D propagators can be easily 
handled, in fact a generic scalar propagator takes the simple form: 

n ^D, a icosxm(n- \y-y'\) ^ ^ , , 

G s (Xrn, y-y) = o : where Xm = vk 2 - m 2 (A. 16) 

2 Xm sin Xmvr 

and m is the 5D mass of the scalar field. One can therefore exploit this form by expanding in 
KK modes along one of the extra dimensions, say Xq, and write the 6D propagator in terms 
of resummed 5D propagators: 

oo 

Gf{h^-^') = Gf( X i,x 5 -x' 5 ) /;(x 6 )/,(A), (A.17) 

l=— oo 

where 



fi{x s ) = -^e iml and X i = Vk 2 - I 2 . (A.18) 
v 2n 



The fi's are the wave functions on a circle and / the KK masses for the 5D modes. This 
method is extremely powerful, especially to calculate corrections for modes like the (n,0): 
the fact that the external fields do not carry any momentum along x 6 , together with the 
orthonormality of the wave functions fi, allows to easily replace the integral in the coordinate 
xq with a sum. The orbifold propagator and the scalar loop we are considering here are given 



by eqs. (A. 4) and (A. 5). Contrary to the winding mode method, this one can be extended 



in a straightforward way to all the other loop diagrams. 
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Torus 

After integrating over x$ and xq, the torus contribution can be written as: 



n T = ^ ± icotXl \ (A.i9) 

l=—oo 

To remove the UV divergence, we decided to renormalized it by regularizing each 5D KK- 
propagators similarly to [llj, so that: 

N f j4 7 1 icotxin — 1 



n T = — / <fk 



1_ £ icotjOTT-l (A 2Q) 



£ = — OO 



After Wick rotation, singling out the I = contribution, one gets: 



Y / ^ poo k% [ coth 

n T = - c(3) + ^ 3 Y, dk 
\ i=i ^° 



7Ti 



2 y/k^+P 



N 

= — A' with A' ~ 1.22; (A.21) 

because of the improper regularization scheme, the finite part is different from the previous 
correct result. However, the structure is the same, thus providing a powerful way to check 
the results obtained with the other methods. 

Glides 

The contribution of the glides is finite, therefore we do not have the same issue with the 
renormalization which arose for the torus one. Under the glide, Xq changes sign and therefore: 

fi{9M) = M-X6 + tt) = (-l) l f-i(x 6 ) ; (A.22) 
due to the orthonormality, the contribution of the / 7^ modes vanishes, and we are left with 

IT G = - f d% C T dx 5 ^o^-\x 5 -9(x 5 )\)sin\nx 5 ) 
4 J Jo 2 xo sin xo* 2tt 2 

After wick rotation, and integrating in £5, we obtain: 

N r°° k 3 N 

n G = -2vr 3 / dk E . E = T 7C(3), (A.24) 
4 Jo kE smhirkE 4 

that agrees with the result obtained using the winding modes. 
Under the second glide 

fi(gr{x 6 )) = Mx 6 + tt) = {-l) l fi(x 6 ) ! (A.25) 
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therefore all modes contribute, and a (—1)' factor appears in the sum: 



AT c r 2 ' K 00 

n G , = - / d 4 k / dx 5 V (-1) 

4 ■/ ^0 



where 



N 

T 



i=— 00 
00 



/■OO 

/ cfe 2tt 2 A;| V (-1)' 

Jo l=-oc 



2fc| 



jicosxK 71 " - \ x 5 - gr(x 5 )\) sin 2 (nx 5 ) 
2 xi sinx«vr 27r 2 

2(fc| + / 2 )+n 2 
(k 2 E + P) (k% + l 2 + n 2 ) 

2kl 



Sinh 7t/c B v / n 2 + fc| sinh ^ v / n 2 + fc 2 ^ 



AT 



[7C(3) + 5 1 (n) 



Si(n) = 2tt 3 / rfp 



Jo 



y/n 2 + p 2 sinh ^Fa/h 2 + p 2 j 



From this expression 



(A.26) 



(A.27) 



Rotation 

Similarly to the glide, for the rotation, only the zero mode contributes: 

IT R = * [d*k f 2 \ X5 ^os Xo (n-\x 5 -r(x 5 )\)s^(nx 5 ) 



4 J Jo 2 xo sinxoTT 2vr 2 

A" f , in 2 

= i,\* k ww^- (A - 28) 

After Wick rotation, we regularize the integral by cut off as before: 

N f A „ 2 , n 2 AT 2 2 , A 2 + n 2 , k . 

R = Jj dk ° 2 * kE WTW) = T n V log — • (A ' 29) 



A. 3 6D Kaluza Klein expansions method 

This method, the most commonly used one, makes use of the expansion in 4D KK modes, 
therefore one needs to compute loops with usual 4D propagators and then sum over the KK 
momenta of the towers. However, computing all the necessary couplings between modes can 
be tedious, and a Fourier transform that goes back to winding modes is necessary for the 
renormalization of the torus contribution. Nevertheless, this method can be easily applied 
to any loop structure. 

Here we will again stick to one concrete example. The contribution of a scalar field with 
parities (p r ,p g ) can be written as 

U Pg Pr = n T + PgU-G + PgPr^-G' + Pr^R 5 (A.30) 
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therefore if we calculated the contribution of all 4 parity possibility, we would be able to 
extract each term: 



n T = 


- (n ++ + n + _ 


+ n_ 4 


_ + n__) 


(A.31) 


n G = 


\ (n ++ - n + „ 


+ n_ 4 


--n__) 


(A.32) 


n<y = 


\ (n ++ - n + _ 


-n^ 


- + n__) 


(A.33) 


n R = 


\ (n ++ + n + _ 


-n^ 


--n__) 


(A.34) 



The couplings that enter the loop, ^(„o)^( n o)^(m()^W are P ro P°rtional to ig 2 rj m with a 
coefficient that depends on the wave function integrals. We listed such coefficients in the 
following table (here m, I ^ are intended): 



(Pr Pg) 


(++) 


(+-) 


(-+) 


(— ) 


(0,0) 


2 








(m, 0) m even 


2 




L 2 




, m^n 
(m, 0) ' 

m odd 




2 




2 


(n,0) 




1 




3 


(0, /) 1 even 


2 






2 


(0, /) 1 odd 




2 


2 




(m, I) m 7^ n 


2 


2 


2 


2 


(n, /) 1 even 


3 


1 


1 


3 


(n, Z) 1 odd 


1 


3 


3 


1 



For example, using the previous table, the correction coming from ++ and 4>+ 
into the loop are: 

iU ++ = N^ 2G(m,Z) +^((-l)'G(n,Z) + 2G(0, 2/)) 

(m,0>l i>l 

+ ^ 2G(2m,0) + 2G(0,0)) , 

m>l,m^ra 

in + _ = iv( 2G(m,Z) +^(-(-l)'G(n,Z) + 2G(0,2Z- 1)) 

(m,0>l i>l 

+ ^ 2G(2m - 1, 0) + G(n, 0)) ; 

where 

/■°° k 3 
GM = J dkE k*+m* + l> 

is the integral appearing in the 4D loop. 
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_ running 

(A.35) 

(A.36) 
(A.37) 



Torus 



For the torus contribution, we reconstruct a sum over all the KK modes on a torus compact- 
ification: following the usual Fourier expansion in the double sum and removing the zero 
winding mode contribution 

^t = N X - G(m,l) = jT 6 . (A.38) 

(m,()ez 2 

Glides 

Following the same procedure, the glide contribution can be written as 

n G = f 2 ^(-l) m G(m,0) = f 7C(3) ; (A.39) 



4 



while for the second glide 



Rotation 

For the rotation contribution, the loop calculation gives: 



N N f A ( 

n fl =-2(G(0,0)-G(n,0)) = T J dk E fk E 



B Appendix: loop integrals 



1.3 



k 2 E + n 2 / 

N , 9l A 2 + n 2 
= -wVlog 2 . A.41 
4 rr 



The functions of n appearing in the loop corrections can be expressed in terms of the three 
following integrals: 

foe l3 

$i(ra) = 2tt 3 / dk (B.l) 

Jo V ^ 2 + n 2 sinh 7r v /c 2 + n 2 

$ 2 ( n ) = 2tt 3 / rfA; v ; , (B.2) 

Jo yk 2 + n 2 sinh7rvA; 2 + n 2 



roc 



fc 3 (VF+?-n) 



$ 3 (n) = 2vr d / dA; — v ' . (B.3) 

n v /c 2 + n 2 sinh 7r v /c 2 + n 2 



34 



Those integrals can be computed analytically, and we found 



$i(n) 
$ 2 (n) 



2im 4> 2 (n) - 03 (n) 
n(7r 3 + -K(j) 2 {n)) , 
1 -2n 2 



2n 



7r d -2 3 (n) + — 4 (n); 



with 



s (n) = Li s (e 



2-7m\ 



- 2 s LLfe 7 ™) - i'- 



V s] 



-n 



s-l 



(B.4) 
(B.5) 

(B.6) 



(B.7) 



where Li s is the Polylogarithmic function of order s and the imaginary term cancels the 
imaginary parts of the Polylogs to make a real function of n. Numerically the integrals are 
suppressed for large n: 





n = 1 


n = 2 


n = 3 


$1 


1.43 


0.109 


0.0067 




0.54 


0.047 


0.0030 


$ 3 


1.02 


0.037 


0.0015 



The functions appearing in the loop corrections are (where we report the numerical value 
for n = 1): 



B 2 
B 3 
Fx 

F-2 



$i ~ 1.43, 

$i - $ 2 ~ 0.89 



9 



-$5 



3$ 3 

$i - 2$ 2 ~ 0.35, 
$ 2 ~ 0.54 . 



1.83. 



(B.8) 
(B.9) 

(B.10) 

(B.ll) 
(B.12) 
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